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Abstract 

Necessary conditions for generating phase controllable asymmetry in spatially symmetric systems 
using lasers are identified and are shown to be identical in quantum and classical mechanics. First, 
by studying the exact dynamics of harmonic systems in the presence of an arbitrary radiation field, 
it is demonstrated that anharmonicities in the system's potential are a necessary requirement for 
phase controllability. Then, by analyzing the space-time symmetries of the laser-driven Liouville 
dynamics for classical and quantum systems, a common set of temporal symmetries for the driving 
field that need to be violated to induce transport are identified. The conditions apply to continuous 
wave lasers and to symmetry breaking effects that do not rely on the control of the absolute phase 
of the field. Known examples of laser fields that can induce transport in symmetric systems are 
seen to be particular cases of these symmetry constraints. 
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I. INTRODUCTION 



Recent years have witnessed the birth and rapid development of the coherent control field 



a, |6j , in which the coherence properties of applied laser fields are employed to steer 
a given quantum dynamical process in a desired direction. Of the different control schemes 
that have so far been developed, there is a general class that has the ability to induce phase 
controllable transport in spatially symmetric systems without introducing a bias voltage in 
the potential, a phenomenon that is referred to as laser-induced symmetry breaking. 



ds 
ll- 



This symmetry breaking effect is typically achieved by driving the system with AC fie^ 
with frequency components nuj and muj, where n and m are integers of different parity 
The nonlinear response of the system to such fields results in net dipoles or currents whose 
magnitude and sign can be manipulated by varying the relative phase between the two 
frequency components of the radiation J?]. For the popular case of n = 1 and m = 2 the 
rectification effect first appears in the third order response of the system to the incident radi- 
ation. At this order the system mixes the frequencies and harmonics of the field, generating 
a phase-controllable zero-harmonic (DC) component in the response. 

Laser-induced symmetry breaking has been demonstrated in a wide variety of systems 
ranging from atoms to solid state samples. Experimentally it has been implemented for 
generating anisotropy in atomic photoionization [8|], symmetry breaking effects in molec- 



ular photodissociation [91 (^^^ ^^^^ ^Qi. lid ), photocurrents in quantum wells [ll] and in 



trinsic semiconductors [l2j, as well as directed diffusion in symmetric optical lattices 13|. 
Theoretically, it has been studied for generating transport in doped [ij] and bulk semi- 
conductors through interband jisl and intraband |l^ excitations, in graphene and carbon 



nanotubes |17| , and molecular wires 



19|, among others. The scenario is of interest since. 



with current laser technology, it can be employed to generate transport on a femtosecond 
timescale. 
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An interesting feature of this laser control scenario is that it has both a quantum H 



201] and a classical 



21 



22| manifestation. Further, the two versions of the effect 
correspond to the same physical phenomenon arising from the nonlinear response of 
material systems to symmetry breaking radiation fields. In this contribution, we isolate 
minimum conditions on the driving field and on the system that is being driven that are 
necessary for the symmetry breaking effect to occur in quantum and classical mechanics. As 
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shown, the minimum requirements in both cases are identical and, further, the effect can 
be accounted for in both mechanics through a single symmetry analysis of the equations of 
motion. 

Specifically, in Sec. [IT] we demonstrate that laser rectification can only occur in systems 
with anharmonic potentials. Subsequently, in Sec. IIIII temporal symmetries of the field that 
need to be violated to induce transport are isolated. This is done by studying the space- 
time symmetries of the Liouville equations of motion for laser-driven quantum and classical 
systems, and by isolating symmetries of the field that rule out any nonzero average currents 



or dipoles at asymptotic times. In doing so we considerably extend a previous analysis [21 1 
that identified conditions on the field necessary for laser-rectification in classical ergodic 
systems. It applies to both quantum and classical systems and makes no assumption of 
ergodicity. 

II. CONDITIONS ON THE SYSTEM 

Consider first the exact solution for the dynamics of a harmonic oscillator in the presence 
of an arbitrary space-homogeneous radiation field E[t). The Hamiltonian of the system 
reads 

H{x,p) = — + -mujQX^ - qE{t)x, (1) 

where x and p are the position and momentum of the particle of mass m and charge q, and 
Uq is the natural frequency of the oscillator. Symmetry breaking here would correspond to 
the production of a net dipole moment. In the quantum case, the dynamics of the position 
xh(^) and momentum operators Ph(^) in Heisenberg picture is dictated by the Heisenberg 
equations of motion 

= \[xnit),Huit)] = -PhW, (2a) 
at ih m 

= l[pH(t), Mt)] = -mco'o^nit) + qE{t), (2b) 

where is the Hamiltonian operator in Heisenberg picture and [/, Hu] = fHn — Hnf for 
any operator /. In the classical case, the position x{t) and momentum p(t) variables obey 
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Hamilton's equations 

^ = {x{t),H{t)} = ^p{t), (3a) 
^ = {pit), Hit)} = -mw'.xit) + qEit), (3b) 

where {f,H} = ^^-^^is the Poisson bracket. The difference between Eqs. (E]) and 
is that the former is a differential equation for operators, with operator initial conditions 
XuiO) = X and Ph(0) = p, while the latter is an equation for functions. 

These two sets of equations can be solved exactly using Laplace transforms. In fact, for 
a general external field of the form 



Eit) = I du eiLu)e'^\ (4) 

oo 



the usual procedure [23| yields: 



xnit) = xniO) cos(ct;ot) + ^ sin(c<Jot) 



, ge(a;) ztu sin(co'ot) + a;o cos(c<Jot) — cuqC*'^* , , 
+ / duj ; (5) 



p(0) 

x{t) = x(0) cos(c<Jot) H sin(u;o^) 

+ n ^^^^-^ sin(u;ot) + c^q cos(^ot) - cupe"^* 

The first two terms in Eqs. (I5l) and (E]) describe the field-free evolution of the oscillator, 
while the third one characterizes the influence of Eit) on the dynamics. 

Note that a driven harmonic system can only oscillate at its natural frequency ujq and at 
the frequency of the field uo. That is, there are no frequency components of the dipole that 
oscillate at multiples or combinations of the frequencies of the field. Hence, if Eit) is unbiased 
(e(0) = 0) no net dipole can be induced. Thus, we conclude that a necessary requirement 
for symmetry breaking in quantum and classical mechanics is that the potential of the 
system is anharmonic. As seen below, the anharmonicities permit the nonlinear response of 
the system to the incident radiation that mixes the frequencies and harmonics of the field 
and, for a certain class of radiation sources isolated below, can lead to the generation of a 
phase-controllable zero-harmonic (dc) component in the response. 
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III. CONDITIONS ON THE FIELD 



We now isolate those temporal symmetries of the field that need to be violated to induce 
transport in both quantum and classical mechanics. To do so we consider a sjTiimetric one- 
dimensional system composed of charged particles coupled to an external field E{t) in 
the dipole approximation. This is done without loss of generality since the polarization of 
the field effectively defines an axis along which symmetry breaking can arise. The system's 
Hamiltonian is then: 

N 2 N 

j=i J j=i 

where and qj denote the coordinate, momenta, mass and charge of the j-th particle 

and X = {xi,X2, ■ ■ ■ ,X]\r). The systems of interest have a potential V{x) that is invariant 
under coordinate inversion V{—x) = V{x), and the driving field E{t + a^) is an arbitrary 
time-periodic zero-mean function, with period T and global phase a. 

In order to keep a close analogy between the quantum and classical case we frame this 



analysis in phase space and adopt the Wigner representation of quantum mechanics [241. 125 1. 
In it, the state of the quantum system is described by the Wigner distribution function 
p^^{x,p,t), which constitutes a map of the system density matrix p in the phase space of 
position X and momentum p variables. For iV-particle one- dimensional systems it is defined 



by |25| 



^) / ■■■/ du^■■■duNeT^P■^{x-u/2\p{t)\x + u/2), (8) 

where \x) = |xi)|x2) ■ ■ ■ \xn) and p ■ u = J^iLiPi^i- In this way the quantum or classical 
Liouville evolution can be expressed as 

Vi3pp{x,p,t) =0, (9a) 

where the label /? indicates either classical (/? = c) or quantum (/3 = W), with pc{x,p,t) 
denoting the classical phase space density. For the Hamiltonian in Eq. (^^, the operator T>/3 



determining the dynamics is given by 

N 



25 



26| 



(9b) 



/ (^^/2)-----^----c^--^--^--n.) 0---^^ 

\i + -+\n=3,5,... 



where the last summation in runs over all positive integer values of Ai, . . . , Aat for which 
the sum Ai + A2 H — ■ + Xn is odd and greater than one. In phase space the formal structure 



of the quantum and classical evolution is remarkably similar |27l. l28l]. In the limit /i ^ the 
second term in Eq. fl^ vanishes and the quantum equation of motion reduces to the classical 
evolution {T>^ — > X'c)- Note that Eqs. (jS]) and ^ are fully consistent with the Hamiltonian 
in Eq. ([7]). However, when the radiation- matter interaction term in the Hamiltonian goes 
beyond the dip ole approximation both of them need to be modified in order to ensure gauge 
invariance 26 1. 



In the absence of an external field the equations of motion [Eq. ([9])] are invariant under 
reflection [x —x, p —p)- Hence, if the system is initially prepared with a given 
phase-space symmetry this initial symmetry is preserved at all times during the subsequent 
dynamics. Symmetry breaking is achieved by coupling the system to E{t). However, if 
E{t) has a zero temporal mean (AC field) then not every E{t) will generate transport. As 
shown below, by properly lowering the temporal symmetry of E{t) it is possible to induce 
rectification in the response. Furthermore, the resulting symmetry constraints on E{t) are 
identical for the classical and quantum case. As will become evident, this is a consequence of 
the important fact that the quantum corrections in "D-yv have the same symmetry properties 
as ^^^TT- under inversion of position and momentum coordinates. 

axj opj 

We focus on rectification effects that survive time averaging and that are independent of 
the global phase a of the laser beam. Typically, symmetry breaking effects that depend on 
a are difficult to control (although not impossible 29|]) since this requires an experimental 
setup that both locks the absolute phase of the laser and has control over the center of mass 
motion with respect to the laboratory frame. Hence, the quantities of interest are the mean 
position and momentum averaged over time and over a: 

dt r-da. 



r' dt da^ / . ^^ . X 

{x)f3=hm / — / — Tr(xp^(a;,p,t)); (10a) 

/■^/^ dt f'^'' da 

{p)^=\im - —TTipp^ix,p,t)y, (10b) 

^^°°J-t/2 t Jo ^tt 

where the double overbar indicates this kind of averaging. Here the notation (■ ■ ■ )/3 denotes 
the classical ensemble average (/3 = c) or quantum expectation value (/3 = W), and the trace 
is an integration over the 2A^-dimensional phase-space {x,p). When the symmetry of the 
system is not broken, both {x)i3 and are zero. Below we determine symmetries of the 
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field and of tlie initial condition that guarantee that this is indeed the case. When these 
symmetries are violated a net dipole or current is expected to appear. 

The fact that we are only interested in a-independent properties eliminates the necessity 
to invoke ergodicity in the analysis. The average over a is sufficient to obviate any initial- 
time preparation effects, which is the main role of the ergodicity assumption in the purely 



classical analysis of Ref. 21 1. 



We now tabulate the symmetries of the field that will be relevant for our purposes. The 
field may change sign every half a period T, 

E(t + T/2) =-E(t); (11a) 

or be symmetric or antisymmetric with respect to some time t' 

E{t - t') = + E{-{t - t'))- (lib) 
E{t-t')=-E{-{t-t')). (11c) 

Each of the symmetries in Eq. f|TT]) leads to a transformation that leaves the equations of 
motion invariant while changing the sign of the position and/or momentum variables. They 
are identical for the quantum and classical case. For example, if E{t) satisfies Eq. ( lllap . 
then is invariant under Ti defined as: 

Ti : t^t + T/2; x ^ -x; -p; (12a) 

where we have taken into account that under inversion of position and momenta, 
X77 A77 and -^^^r- X77 ^ X77 smce Ai + A2 + ■ — h Aat m Eq. 09d 

dxl^-dx)l' dx\^-dx]f dpl^-dp%'^ dp\^-ap]f ^ ^ ^ 

is odd. Similarly, if E{t) satisfies Eq. (lllbl) [or Eq. flllcp ] then Vjj is invariant under T2 [or 
73], where 

T2 : t-t' ^ -it-t')] x^x; p^-p; (12b) 
%: t-t' ^ -it-t'); x^-x- p^p. (12c) 

Other temporal symmetries of the field exist but play no role in this analysis since they do not 
lead to invariance transformations that change the sign of the position and/or momentum 
variables. 

Now, given a solution to Eq. ([9]), p[^{x,p,t), if Vj3 is invariant under %, one can generate 
another solution to the same equation by applying 7^ to pp{x,p,t). The new solutions 



Pp {x,p,t) = TaPf3{x,p,t) generated by the invariance transformations in Eq. f|T2l) are: 

pJJ^(a;,p,t) = Tipp{x,p,t)= pi3{-x,-p,t + T/2); (13a) 
pf{x,p,t)= T2Pf3{x,p,t -t') = pp{x,-p,~{t-t'))- (13b) 
pyix,p,t)= %pf3{x,p,t-t') = pp{-x,p,-{t-t')). (13c) 

Further, if the original solution pp{x,p,t) predicts an average position {x)/3 and momenta 
(p)/3, the transformed solutions p^^\x,p,t) will predict a mean position (^)^"^ and/or mo- 
menta {p)^0^ that has the same magnitude but is opposite in sign: 



= -(2;)/3; 


(P)? = 


-(P)/3; 


(14a) 




{P)f = 


-(P)/3; 


(14b) 




= 




(14c) 



The argument is completed by showing that the average position and momenta predicted 
by Pj3{x,p,t) and p'^^\x,p,t) coincide. If this is the case, it follows from Eq. f|T^ that 
symmetry breaking cannot occur. For this we exploit the possible symmetries of the initial 
state: 

P/3(a?,P,to) = P/3(-a;, -p,^o); (15a) 
Pp{x,p,to) = pp{x,-p,to)] (15b) 
pp{x,p,tQ) = pp{-x,p,tQ). (15c) 

The first one corresponds to a state with zero mean position and momenta, while the second 
and third describe an initial state with either zero mean momenta or zero mean position, 
respectively. 

Consider the case in which the equations of motion are Ti invariant. The distributions 
Pi3{x,p,t) and py{x,p,t) satisfy the same equation of motion but do not, in general, coin- 
cide. However, if the initial condition for the original solution p^(x,p, to) is invariant under 
reflection [Eq. f ll5al) ] then 

Pp\x,P,to-T/2) = pp{-x,-p,to) = pfi{x,p,to) = py{-x,-p,to-T/2). (16) 

That is, the original and transformed solutions start from the same initial distribution but at 
initial time they experience a different value for the global phase of the field, E{to + a^) and 
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E{to + {a—7i)j^) = —E{to+a^) respectively. Since the averages in Eq. ffTOl) are independent 
of a, they coincide for the two solutions. Hence, no rectification can be induced when the 
field satisfies Eq. f ilial) and the initial condition satisfies Eq. fllSal) . 

The argument for the two other cases is very similar. If the field satisfies Eq. ( Illbj) [or 
Eq. (11 Id) ], the equation of motion is 7^ (or Tj) invariant. Even when the original pp{x,p,t) 
and transformed p^p\x,p,t) [or p^^\x,p,t)] solutions obey the same equation of motion, 
they do not need to coincide. However, if the initial condition of the original solution satisfies 
the symmetry in Eq. fll5b|) [or Eq. (I15cl) ]. then 



pf{x,p,-to + t') = pp{x,-p,to) = pp{x,p,to) = pf{x,-p,-to + t') (17) 
pf{x,P,-to + t') = pp{-x,p,to) = pp{x,p,to) = pf{-x,p,-to + t') (18) 

The transformed solution has the same initial condition as the original one but as we had 
prepared the system a time 2to — t' before. The difference between the two solutions is that 
they experience a different global laser phase at preparation time. Since we are not interested 
in effects that depend on the global laser phase, the mean position and momenta in Eq. (iTOll 
for the original and transformed solution need to coincide. However, from Eq. ( 114b[) [or 
Eq. (I14cl) ] we conclude that this can only happen if (p)/3 = [or {x)i3 = 0]. 

In summary, for spatially symmetric classical or quantum systems initially prepared in 
a symmetric state that satisfies Eq. f|T5l) . net transport using time-periodic external fields 
with zero temporal mean can only be generated if the field violates the temporal symmetries 
in Eq. ( ITTi) . Further, any symmetry breaking effect that may be achieved with a field that 
satisfies Eg. (ITTl) is necessarily due to an effect that depends on the global phase of the laser 



(cf. Ref. |29(). 

It is natural to ask what kind of fields have sufficiently low temporal symmetry to induce 
net transport. Monochromatic sources satisfy all the symmetries in Eq. (ITT]) and, as a 
consequence, cannot be used to induce symmetry breaking. However, by adding a second 
frequency component to a monochromatic source it is possible to lower the symmetry of the 
field and induce symmetry breaking. For instance, a field like 

E(t) = enuj cos{nujt + 0„^) + e^^ cos{mujt + (j)^^), (19) 

where n and m are coprime integers so that E{t) has a period T = 2tt/uj, satisfies Eq. (fTTj) 
only under special conditions. It satisfies (I flap only if n and m are odd. Thus, a field with 



n = 3 and m = 1, like the one used in the 1 vs. 3 photon control scenario [ij, will not be 
symmetry breaking. However, a field with n = 2 and m = 1, like the one employed in the 
1 vs. 2 scenario, does not satisfy Eq. (ITTl) and is expected to induce net dipoles and cur- 
rents. These dipoles and currents are phase-controllable since, by varying the relative phase 
between the two components of the beam, the uo + 2uj field may satisfy Eq. (]llb|) or flllcp 
and thus rule out the possibility of inducing currents or dipoles, respectively. Explicitly, 
when (j)2u) — 20a; = 0, ±7r, ±27r, ■ ■ ■ an u; + 2ti; field satisfies Eq. flllbl) and zero currents are 
expected. Similarly, when 02^ — 200; = ±f , ■ ■ ■ it fulfills symmetry (lllcp and no dipoles 
can be induced. For all other choices of the phases an a; + 20; field is expected to induce 
symmetry breaking. 



IV. CONCLUSIONS 



In conclusion, we have shown that the minimum conditions for the generation of phase 
controllable asymmetry in spatially symmetric quantum and classical systems using time- 
periodic external fields with zero temporal mean are identical: anharmonicities in the sys- 
tem's potential are required as is a driving field that violates the temporal symmetries in 
Eq. (fTTj) . These conditions refer to symmetry breaking effects that do not rely on the control 
of the absolute phase of the field. The derived results provide necessary conditions for the 
generation of asymmetry, applicable to all systems. Additional sufficient conditions may be 
required, but these depend upon the specific system under consideration. 

Further, we have shown that both quantum and classical versions of the symmetry break- 
ing effect can be accounted for through a single space-time symmetry analysis of the equa- 
tions of motion. The anharmonicities in the potential permit the nonlinear response of 
the system to the incident radiation that, through harmonic mixing, and for fields that 
violate Eq. flTT]) . can lead to the generation of a phase-controllable DC component in the 
photoinduced dipoles or currents. 
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